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表面 活性 剂 溶液 粗 粒 化 模型 couette 流动 模拟 
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摘要 : 本文 主 要 利用 粗 粒 化 分 子 动力 学 方法 对 不 同 浓度 下 表面 活性 剂 水 溶液 在 couette 流动 过 程 中 胶 束 的 形成 、 速 度 分 布 与 流 
变性 规律 进行 了 研究 。 结 果 表明 ， 对 于 相同 浓度 的 表面 活性 剂 溶液 ， 随 着 剪 切 速率 的 增 大 ， 胶 束 完全 动态 平衡 时 间 越 来 越 
小 ; 而 随 着 表面 活性 剂 分 子 请 度 的 增加 ， 胶 束 动态 平衡 时 间 也 越 来 越 小 ， 形 成 胶 束 越 来 越 快 。 在 剪 切 速率 与 溶液 剪 切 粘度 关 
系 方面 ， 保 持 表面 活性 剂 分 子 浓度 不 变 情 况 下 ， 随 着 剪 切 速率 的 增 大 ， 剪 切 粘度 逐渐 减 小 ， 呈 现 剪 切 稀 化 特征 。 
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A coarse-grain molecular dynamics simulation for surfactant aqueous solution in 


couette flow 
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Abstract: In this paper, the formation of micelle, velocity distribution, and rheological behavior in surfactant 
aqueous solution under couette flow were studied using coarse-grained molecular dynamics simulations. Simulation 
results show that the dynamic equilibration time becomes shorter with the increase of shear rate on the condition of 
the same surfactant aqueous solution concentration. If the shear rate is a constant number, the formation of a micelle 
was more rapid when the surfactant molecular concentration increases in the aqueous solution. When the surfactant 
molecular concentration keeps constant, the shear viscosity becomes smaller with the increase of shear rate, and the 
shear-thinning characteristic can be obtained. 
Key words: surfactant aqueous solution; Coarse-grained molecular dynamics simulations; couette flow; shear 
viscosity 


0 引言 


随 着 我 们 对 能 源 的 需求 量 越 来 越 大 ， 新 能 源 的 学 模拟 研究 ， 但 是 具体 到 双亲 性 表面 活性 剂 ， 在 表 
开发 以 及 能 源 的 高 效 利用 就 变 得 越 来 越 重要 ， 而 减 面 活性 剂 溶液 分 子 动力 学 模拟 方面 的 研究 较 少 ， 
阻 作为 一 种 有 效 减 小 管道 内 流动 阻力 的 方法 ， 可 以 N.A. Volkov, N. Arai 等 人 的 研究 比较 有 代表 性 。 其 


达到 节能 降 耗 的 目的 。 所 以 减 阻 用 表面 活性 剂 的 研 中 N.A. Volkov[8] 对 十 二 烷 基 硫酸 钠 (Sodium 
完成 为 目前 的 研究 热点 之 一 04。 而 表面 活性 剂 在 溶 dodecyl sulphate, SDS) 溶液 中 胶 束 形成 进行 分 子 动 
液 中 的 微观 结构 是 影响 表面 活性 剂 流 变性 的 关键 ， 力学 模拟 ， 发 现 NaCl 的 加 入 对 形成 大 尺寸 胶东 更 
因此 深入 了 解 表面 活性 剂 在 溶液 中 的 胶东 形成 和 流 有 利 ， 该 模拟 详细 描述 了 胶东 的 形状 特性 。N. 
变性 质 具有 重要 意义 。 Arai9] 采 用 DPD 方法 研究 了 纳米 尺度 平行 通道 内 表 
目前 对 于 表面 活性 剂 溶液 的 研究 主要 集中 在 实验 测 面 活性 剂 溶液 在 剪 切 流 场 中 的 流 变 特性 ， 发 现 通过 
量 阶段 ， 即 通过 宏观 流 变 性 的 测量 来 推测 微观 结构 对 纳米 平面 通道 内 表面 进行 疏水 - 亲 水 化 学 改 性 ， 可 


或 者 直接 通过 电镜 对 表面 活性 剂 的 微观 结构 进行 观 以 使 表面 活性 剂 洲 液 实现 由 剪 切 变 稀 - 剪 切 增 稠 的 改 
察 67。 在 微观 结构 模拟 方面 ， 一 般 都 是 针对 普通 聚 变 。 而 A. SambasivamIt19 则 首次 采用 分 子 动力 学 
合 物 的 胶 束 形成 、 断 裂 、 流 变性 等 方面 的 分 子 动力 (Coarse-grained molecular dynamics, CGMD) 模拟 
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方法 对 CTAC+NaSal 溶液 剪 切 流 场 下 的 棒状 胶 束 变 
形 和 断裂 行为 进行 模拟 ， 通 过 能 量 分 析 方 法 对 单 体 
胶 束 长 度 和 剪 切 流 场 下 胶 束 取向 分 布 进行 了 研究 ， 
其 中 得 出 了 CTAC 胶 束 长 度 平均 值 为 5.4nm 的 结 
论 ， 该 模拟 的 总 时 间 长 度 为 3 bs。 在 胶 束 形成 时 间 
和 胶 束 大 小 的 定量 研究 方面 ，Wataru “Shinodall0 模 
拟 了 聚 氧 乙烯 等 几 种 表面 活性 剂 在 空气 -水 、 油 -水 
界面 上 胶 束 形成 过 程 ， 对 胶 束 大 小 与 胶 束 形成 时 间 
进行 了 定量 计算 。 在 粗 粒 化 表面 活性 剂 溶液 couette 
流动 在 大 剪 切 率 情况 下 流 变 性 研究 方面 ，J Castillo- 
Tejas023 对 十 六 烷 基 三 甲 基 氧 化 铵 (CTAC) 表面 活 
性 剂 溶液 的 couette 流动 和 poiseuille 流动 进行 了 粗 
粒 化 非 平衡 分 子 动力 学 模拟 ， 并 且 首 次 对 CTAC 单 
体 进行 了 一 个 杀 水 头 基 和 四 个 跪 水 尾 基 的 粗 粒 化 假 
定 ， 得 到 了 不 同 表面 活性 剂 浓度 下 溶液 couette 流动 
剪 切 速度 与 剪 切 应 力 的 关系 ， 剪 切 速度 与 剪 切 粘度 
的 定量 关系 。 作 为 对 J。，Castillo-Tejas 模拟 应 用 范围 
WE, E.O. Castrejbn-Gonzdlez03] 实 现 了 水 -表面 
活性 剂 分 子 - 油 溶液 混合 物 在 大 剪 切 率 情况 下 couette 
流动 的 粗 粒 化 分 子 动力 学 模拟 ， 得 到 了 随 着 剪 切 速 
率 的 增加 ， 洲 液 呈 现 由 牛顿 流体 向 剪 切 变 稀 ， 再 向 
剪 切 增 稠 ， 甚 至 最 后 再 向 牛顿 流体 转变 的 结论 。 

综 上 所 述 ， 在 胶 束 形成 时 间 的 研究 方面 ， 研 究 
者 通常 使 用 比较 精细 的 分 子 模型 进行 胶 束 形成 研 
究 ， 但 是 这 种 研究 有 两 个 缺点 ， 即 通用 性 差 和 计算 
量 巨大 。 所 以 本 文采 用 粗 粒 化 分 子 动力 学 方法 进行 
研究 ， 简 化 了 分 子 模拟 所 考虑 的 分 子 极 性 和 不 同 基 
团 之 间 的 受 力 ， 得 到 胶 束 的 形成 时 间 。 在 对 表面 活 
性 剂 溶液 的 研究 方面 ， 与 前 人 的 研究 相 比 ， 主 要 多 
考虑 了 溶液 浓度 和 温度 这 两 个 变量 ， 对 表面 活性 剂 
的 流 变 性 进行 了 定性 分 析 ， 得 出 了 与 目前 宏观 实验 
结果 相近 的 结论 ， 加 深 了 对 表面 活性 剂 溶 液 流 变性 
的 认识 。 
1 数值 方法 

本 文中 坐标 xyz 方向 的 表面 活性 剂 溶液 模拟 区 
域 分 别 为 300 X200 X150 (上 下 壁面 区 域 Y 坐 
标 分 别 为 20c ~22.50 和 -2.5o ~0o ) ， 其 中 每 个 
粗 粒 化 表面 活性 剂 分 子 由 一 个 亲 水 头 基 (直径 为 
4o) 和 四 个 朴 水 尾 基 〈 直 径 为 c ) 组 成 ， 粗 粒 化 
水 分 子 直径 也 为 o 。 其 中 表面 活性 剂 分 子 头 基 - 表 面 
活性 剂 头 基 、 表 面 活性 剂 头 基 - 粗 粒 化 水 分 子 、 表 面 
活性 剂 尾 基 - 表 面 活性 剂 尾 基 、 粗 粒 化 水 分 子 - 粗 粒 


化 水 分 子 之 间 的 作用 均 采 用 LJ (12-6) 模型 ， 而 表 
面 活性 剂 头 基 - 其 它 表 面 活性 剂 分 子 尾 基 、 表 面 活性 
剂 尾 基 与 粗 粒 化 水 分 子 之 间 采 用 LJ9-6 模型 施加 排 
斥 力 ， 如 图 1 所 示 。 
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图 1 表面 活性 剂 与 水 分 子 受 力 示意 

Fig.l Surfactant and water CG molecule interaction 
表面 活性 剂 分 子 头 基 与 相 邻 尾 基 、 尾 基 与 相 邻 
尾 基 之 间 采 用 Harmonic 势能 模型 。 表 面 活 性 剂 溶 
液 采 用 NVT 系 综 ， 固 体 壁 面 采 用 NVE 系 综 。 粗 粒 
化 表面 活性 剂 溶液 couette 流动 如 图 2 所 示 。 
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图 2 粗 粒 化 表面 活性 剂 溶液 couette 流动 示意 
Fig.2 The surfactant aqueous solution in couette flow 

1.1 L-J 模型 

L-J (12-6) 势能 表达 式 : 


12 6 
o o 
U; =42 z) {z| +8,%, SY, 
r, To 
y y 


截断 半径 x=2.50 。 
排斥 势能 函数 (Replusive potential) : 


9 
1.050, 
rj 


可 以 用 LJ96 模型 近似 替代 。 
1.2 Harmonic 模型 
Harmonic 势能 函数 表达 式 : 
U, (r) = Ky, =a 

ho j 通过 Lorentz-Berthelot 法 则 : 
o= 0.5(o, + 0 ) 计 算得 到 。 
1.3 NVT 与 NVE 系 综 

NVT ENZ (EH Nose-Hoover 恒温 条 件 ) 
的 方程 为 : 


dr 

dt 

dy 万- 
dt m, ~ 
dé 1 


r, vyo FRAP i WME REM, TA 
目标 温度 ; LAA RE: & 为 热 动力 学 摩擦 系 

A, 而 Q 定义 为 恒温 质量 ， 表 达 式 为 O=LTr?，T 
为 系统 的 松弛 时 间 。 


布 ， 排 布 后 的 各 个 粗 粒 化 分 子 坐 标 信息 作为 模拟 的 


初始 条 件 。 


况 运行 150 万 步 ， 约 为 


此 非 平衡 态 分 子 动力 学 “NEMD) 模拟 
中 ， 各 工 况 的 无 量 纲 时 间 步 长 均 为 0.001， 每 组 工 


5.5nso 


表 1 各 物理 量 无 量 纲 化 公式 


Table 1 Dimensionless of parameters 


物理 量 名称 


术 
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t=t'r=t'(mo"/e)"” 
p=pmo” 
f=T a, 


1/2 


y=y (o’e/m) 


f= f eo" 
p=p eo" 


3 模拟 结果 与 讨论 


3.1 胶 束 形成 时 间 分 析 


通过 对 胶 束 形成 过 程 的 可 视 化 观测 《使 用 Ovito 


开源 软件 ，Open 


visualization 


tool) ， 可 以 得 出 


couette 流动 前 切 速率 分 别 为 0.03、0.06、0.1、 


NVE 系 综 主 要 在 壁面 使 用 ， 系 统 的 原子 数 、 能 


量 、 体 积 保持 不 变 。 
1.4 其 它 模拟 细节 


各 个 粗 粒 化 分 子 使 


| Material Studio 进行 初始 排 


0.15, 0.20. 0.25 六 种 工 况 下 表面 活性 剂 粗 粒 化 分 
子 形成 胶 束 的 时 间 。 胶 束 形 成 的 典型 过 程 一 般 为 初 
始 状 态 -小 胶 束 团 状态 -大 胶 束 团 状态 -完全 动态 平衡 
状态 。 以 表面 活性 剂 分 子 浓度 0.20， 剪 切 速率 0.1 


为 例 进行 说 明 ， 分 别 如 


图 3a，3b，3c，3d 所 示 。 
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(c) 大 胶 束 团 状态 0. 73ns 


Chin 


(d) 完全 动态 平衡 状态 3. 28ns 


图 3 胶 束 动 态 平衡 状态 


Fig.3 The dynamic equilibrium state of micelles 
度 方向 ) 从 0c 开始 到 200 结束 划分 出 100 个 等 距 
离 长 方 体 空间 ， 即 每 个 小 长 方 体 空 间 的 体积 为 30G 
X0.20 X150 。 对 每 个 切片 长 方 体 空间 内 的 粗 粒 
化 粒子 (包括 表面 活性 剂 分 子 和 水 分 子 〉 进 行 一 段 
时 间 内 的 速度 平均 计算 系 综 达 到 动态 平衡 时 各 个 
切片 空间 内 的 平均 速度 。 


表面 活性 剂 浓度 分 别 为 0.05、0.15、0.20 下 不 
同 剪 切 率 与 胶 束 形成 时 间 的 关系 对 比 ， 如 图 4 所 
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图 4 剪 切 速率 与 胶 束 形成 时 间 
Fig.4 shear velocity and micelle formation time 

可 以 看 出 ， 表 面 活性 剂 分 子 浓度 分 别 为 0.05、 
0.15、0.20 情况 下 ， 随 着 couette 流动 上 板 剪 切 速 
率 的 增 大 (前 切 速率 分 别 为 0.03、0.06、0.1、0.15、 
0.20、0.25)， 胶 束 完全 动态 平衡 时 间 越 来 越 小 。 而 
在 相同 剪 切 速率 情况 下 ， 随 着 胶 束 浓度 的 增加 ， 胶 
束 动态 平衡 时 间 也 越 来 越 小 ， 即 胶 束 形成 速度 也 越 
来 越 快 。 
3.2 表面 活性 剂 浓度 与 溶液 速度 分 布 
研究 表面 活性 剂 分 子 形成 的 胶 束 对 couette 流动 
流 场 的 影响 ， 以 couette 流动 前 切 速率 为 0.1、 胶 柬 
浓度 分 别 为 0、0.05、0.15、0.20 的 速度 分 布 为 研究 
对 象 进行 说 明 ， 得 到 胶 束 浓度 与 速度 分 布 关系 ， 如 
图 5 所 示 。 具 体 实现 办 法 : 在 立方 向 上 《〈 即 速度 梯 
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图 5 胶 束 浓度 与 空间 平均 速度 分 布 


Fig.5 The relationship between micelle concentration and the 


spatial mean velocity distribution 


从 图 5 可 以 看 出 ， 在 剪 切 速率 为 0.1 情况 下 ， 


水 溶剂 (表面 活 怕 


E 剂 浓度 为 0， 图 示 中 黑 线 表 示 ) 


的 速度 分 布 基本 是 单调 递增 的 ， 成 斜率 为 1 的 线性 
分 布 。 在 加 入 粗 粒 化 表面 活性 剂 分 子 以 后 ， 在 立方 
向 0-11 o 区 域 平均 速度 大 于 水 溶剂 ， 而 11-200 区 


域 速度 则 小 于 纯 水 。 至 于 是 否 表 面 活 性 剂 浓 度 越 


小 ， 速 度 分布 1 


HAZE KB), BETA 


持 ， 因 为 表面 活性 剂 分 子 浓度 为 0.1 
区 域 与 浓度 为 0.20 时 相 比 更 加 偏离 纯 水 的 速度 分 


FPF 的 数据 不 予 支 
5 时 在 11-200 


布 。 所 以 可 以 得 出 一 个 结论 ， 表 面 活性 剂 浓度 大 小 
与 改变 couette 流动 速度 分 布 的 幅度 大 小 没有 必然 关 
系 。 
3.3 ARR SIME KA 

剪 切 粘度 由 剪 切 应 力 与 剪 切 速率 得 到 ， 计 算 公 
式 为 : 


P 


{== 
7 
其 中 为 剪 切 粘度 ， 忆 ,为 剪 切 应 力 ，7 为 前 
切 速率 ， 本 文 计算 中 都 进行 了 无 量 纲 化 。 

无 量 纲 温度 T=1. 35、 无 量 纲 剪 切 速率 六 分 别 为 
0.01, 0.025. 0.05, 0.1, 0.15, 0.2, 0.25 F, BYY) 
速率 与 剪 切 应 力 、 剪 切 粘度 的 关系 ， 分 别 如 图 6、 
图 7 所 示 。 
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6 剪 切 速率 与 剪 切 应 力 关 系 


Fig.6 The relationship between shear rate and shear stress 
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图 7 剪 切 速率 与 前 切 粘度 关系 


Fig.7 The relationship between shear rate and shear viscosity 

从 图 6 可 以 得 出 如 下 规律 : 随 着 剪 切 速率 的 增 
大 ， 剪 切 应 力 随 之 增加 。 而 随 着 粗 粒 化 表面 活性 剂 
分 子 浓度 增 大 ， 剪 切 应 力也 越 大 。 
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从 图 7 可 以 看 出 ， 除 了 和 剪 切 速率 为 0.01 和 
0. 025 时 剪 切 粘度 分 布 有 些 杂 乱 ， 没 有 呈现 出 剪 切 
稀 化 的 特征 之 外 ， 之 后 随 着 剪 切 速率 的 增 大 ， 剪 切 
粘度 逐渐 减 小 ， 呈 现 剪 切 稀 化 特征 。 剪 切 速 率 太 小 
(<0.03) ， 粗 粒 化 分 子 的 热 运动 干扰 过 大 ， 导 致 
系统 内 各 层 粗 粒 化 分 子 不 能 完全 做 有 序 的 couette 流 
动 ， 所 以 只 取 0. 05-0. 25 有 效 剪 切 速率 区 域 进 行 分 
析 。 同 时 可 以 看 出 粗 粒 化 表面 活性 剂 分 子 浓度 与 剪 
切 粘度 之 间 呈 现 正 相 关 ， 即 表面 活性 剂 浓度 越 大 ， 
剪 切 粘度 越 大 。 


4 结论 
本 文 主要 对 不 同 浓度 的 表面 活性 剂 溶液 couette 


流动 进行 了 非 平 衡 分 子 动力 学 模拟 ， 得 出 couette 流 
动 中 剪 切 速率 与 粗 粒 化 表面 活性 剂 分 子 胶 柬 形成 时 
间 、 表 面 活性 剂 浓 度 与 溶液 速度 分 布 、 剪 切 速率 与 
剪 切 应 力 、 剪 切 速率 与 剪 切 粘 度 的 关系 。 发 现 表 面 
活性 剂 浓度 和 剪 切 速率 与 胶 束 形成 时 间 负 相关 ， 表 
面 活 性 剂 浓 度 与 溶液 速度 分 布 无 必然 联系 。 在 浓度 
一 定 情况 下 ， 剪 切 速率 越 大 ， 剪 切 应 力 越 大 ， 剪 切 
粘度 越 小 ， 呈 现 剪 切 稀 化 特征 。 通 过 对 表面 活性 剂 
溶液 中 这 些 特定 参数 的 分 析 研 究 ， 可 以 更 深入 地 了 
解 表面 活性 剂 溶液 胶 束 的 形成 以 及 胶 束 的 流 变 特 
性 。 
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